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^\ ' We propose a new and general method for deriving exact density functionals in one dimension 

On , for lattice gases with finite-range pairwise interactions. Corresponding continuum functionals are 

derived by applying a proper limiting procedure. The method is based on a generalised Markov 
^ ' property, which allows us to set up a rather transparent scheme that covers all previously known 

5^ , exact functionals for one-dimensional lattice gas or fiuid systems. Implications for a systematic 

■ construction of approximate density functionals in higher dimensions are pointed out. 

o 

(N 

, Density functional theory provides one of the most powerful methods for studying thermodynamic properties of 

O ' inhomogeneous fluids ||l|-^. Various schemes are known to approximate the excess part of free energy functionals, 

^ which contains the contribution from the interactions between the molecules. Exact density functionals, however, have 

y only been derived in one dimension d — 1. These functionals have their merits in establishing "boundary conditions" 

^ in the sense that approximate functionals used in d > 1 should become exact upon dimensional reduction to d = 1. 

Ci In this Letter we will present a simple general method how to derive exact density functional in one dimension for 

c/2 systems of particles interacting via a finite-range pair potential. In particular we will show that density functionals 

^ , calculated previously |p[-p^ can be recovered. We will first consider lattice systems from which the corresponding 

^ ■ continuum functionals can be straightforwardly derived by a proper limiting procedure. 
The lattice Hamiltonian reads 



. H ^^XsUs +^Vij X^XJ , (1) 

O , where Vij is the interaction potential, Us = ~ fj, is the external potential minus the chemical potential /i, and Xs 
are the occupation numbers {xs = 1 if site s is occupied and Xs = else). A particularly simple situation arises for 
[ a system composed of hard rods of length 2to that is described by an interaction potential vtj = oo for j — i < 2m 
J> • and Vi^j — else. For this model an exact density functional was derived by Robledo and Varea ||] (see also [^). 

The corresponding continuum hard rod fluid has been treated by Percus, who succeeded to set up the first exact 
\l functional for this many-particle system with athermal interactions Later exact functionals have been derived 
for the multicomponent hard- rod and sticky core fluid as well as for the Takahashi model ||l^. The perhaps 
^ J most general functional was set up for a lattice system as in (|l]) with an arbitrary interaction potential exhibiting a 
' finite range |l^ (w^.j =0 for j — i>l)- This system, that we will refer to as the "finite-range model" in the following, 
^\ ' includes as a special case the Ising model for which various exact functionals have been derived separately for periodic 
boundary conditions ||l^-[l^ . 

' We will focus here on the finite-range model with fixed boundary conditions, that means we consider a finite 
^ , lattice of M sites (s = 1, . . . , M) and fix the occupation numbers at the boundary points s = and s = M -I- 1 to be 
I xo = xm+1 = 1- Our basic idea for deriving the density functional is to generate the joint probability distribution 
^ , x[^1t ■ ■ t^m\ by a generalised Markov process, and then to equate the resulting expression with the Boltzmann 
g ' formula x\^ii ■ • ■ ; ^m] = exp(— H) (for convenient notation we set k^T — 1). 

^ Let us start by expressing the joint probability distribution in terms of the conditional probabilities 

• • ■ UJs[Xs\Xs-l, ■ ■ ■ ,Xi], 



. x[xi,---,Xm]^ \_YuJs[Xs\Xs-l,---,Xi\. (2) 

Due to the finite range of the interaction, the conditional probabilities fulfil the Markov property lOs[xs\xs-i, ■ ■ ■ , xi] — 
ujs[xs\xs-i, ■ ■ ■ ,Xg_(^i_i^] ( [|l5|; we set Xs = for s<0 for convenient notation). Factorising Ws[xs|xs_i, . . . , 



into 2' terms corresponding to the possible configurations of {xs, ■ ■ ■ ,Xs_n^i\) yields 



i-i 

UJs[Xs\Xs-l,---,Xs-{l-l)]^'[[ n 

'•=0 jeAs, 



. i-xsi n n (1 



(3) 
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ai = E[xs n n (1 - ^^)] ' = nu n (1 - ^^)] ■ 

E[. . .] denotes an expectation value with respect to x ^-i^d As,r is the set of ah configurations (a;5_(;_i), . . . , Xg^i) with 
exactly r sites being occupied. The superscript j specifies that particular element of As,r, for which the r sites ji , . . . , 
are occupied, and /J^ is the corresponding set of all the remaining vacant sites. Note that a^/bj is the probability 
that the site s is occupied under the condition that the r sites ji , . . . , jr are occupied while the remaining sites i £ ^ 
are empty. Although eq. (j^) looks complicated, its structure becomes immediately obvious when considering the 
special case I = 2, where ujs[xs\xs-i] = u;s[l|l]''='==-it^Jl|0]^=(i-^=-i)ws[0|l](i-="=)='=-iwjO|0](i-^=)(i-^=-i' with = 
'E[xsXs~i]/'K[xs^i\ etc. 

A grand canonical functional Q, = — log Z is now derived by inserting eq. (^) into eq. (^) , and by using for x 
Boltzmann formula. After taking the logarithm and the expectation value on both sides of eq. (0), we obtain 



M M l-l 

s=l 0<i<j<M+l s=l r=0 je^s,, 



log -f +(6J-ai)log 



bl 



(4) 



where for r > 1 and ki < . . . < kr we define pk^ kr = '^[xki ■ ■ ■ Xk^]- According to the definitions of aj and 6J (see 

eq. @), is a functional of all moments Pki,...,kr (with k,. — ki < I), that means it can be considered as a functional 
of (local) correlations up to order Z — 1. In order to derive a density functional that depends only on the occupation 
probabilities Ps, the correlations have to be expressed in terms of the Ps- According to the theorem of Mermin 
this can be done in a manner that the corresponding "correlator equations" are invariant with respect to Ug and we 
will show next that there exists a general method to achieve this . 
To illustrate the method, we choose the special interaction potential 

{oo , < j — i < 2m, 
Vij , 2m <J — — 1< 4m . 

This potential describes hard rods with lengths 2m, which interact via the potential Vij, and it corresponds to a 
lattice version of the continuum Takahashi model (for Vij — Vj^i, p^). The special form of Vij implies that in an 
interval of size 1^1 there can at most two occupation numbers be equal to one. Accordingly, af, and bf, depend only 
on the occupation probabilities {ps} and the two-point correlations 

In order to calculate the two-point correlations, we consider the following three cases: (i) all Xs {s — 1,..., M) are 
zero, (a) exactly one of the Xs is equal to one (and the rest equal to zero), and (Hi) exactly two occupation numbers 
Xi, Xj with 2m <j — i<l are equal to one. By expressing the probabilities for these cases in terms of the Boltzmann 
formula, we find (i) xi^-U Xs = 0] = x[0] = (H) x[xs = l;rest 0] = x[xs = 1;0] = 2^""'^ exp(— u^), and (in) 

x[a;i = 1, a;^ = 1; 0] = 2^^exp[— — iij — Vi,j\ p9| ]. On the other hand, by considering eqs. (0J3) for the cases (i)-(iii) 

we obtain expressions containing the sets {ps} and {pi j}. For example, x[0] = Y[sLii^~''^'^ /b^s^^)- The correlator 
equations then follow by multiplying the expressions for the cases (i)-(iii) in such a way that Z and {u^} cancel, i.e. 
in the present situation we get x[0]x[^i — Ij^^j = 1;0] — ^'^'P{~^i,j)x[xi — l;0]x[^j — 1;0]- Inserting the x's from 
eqs. (HI) gives, for i -I- 2m<j<i + 1 — 1, 



j—l j~2m 

= E[xjY[i^-Xk)] =Pj-Yl Pi^^ ' = Pi^o ' (6) 



k=]-{l-l) k='j-(l-l) 



j — l J 1 J — 1 — 2 m j — l J ~ 1 

bf ^nW {^~^k)] = i-Y. Pk+ E pf^'-, f^^pr- E p^^^- 

fc=j-(i-l) fe=j-(i-l) k=j-{l-l) r=k+2m k=i+2m 

Note that the method for determining the two-point correlator equations can also be applied for deriving higher order 
correlator equations. These equations would be needed for more general interaction potentials. 

Having derived the correlator equations, we can view the moments in the functional fl (^ as functions of the Ps 
and minimise f2 with respect to the Ps to obtain the "structure equations". In fact, these equations can be derived 
more directly by using the relation x[0] exp(— Ui) — x[xi — 1; 0], which gives 
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(Equation (^ together with the correlator equations (^) determines the occupation probabihties Ps-) 

An exphcit solution of the correlator equations (|^) can be obtained for / — 2m + 1 , which means that two rods only 
interact if they are in contact (we set Vi = Vi^i+2m then). The solution is 



Pi,i+2m = 



Pk) +Pi +Pi+2r. 



2[e*' - 1] 



4[e''* - l]piPi+27n 



Pk)+Pi+Pi+2m 



(8) 



After inserting (^ into (^ one also obtains the structure equations in explicit form. For to = 1/2 and Vi — ^Ji in 
particular, the system is equivalent to the Ising model and the structure equations simplify to the ones previously 
derived [p0|-p2[ by using very different approaches. 

Unfortunately, for / > (2m + 1) we are not able to give an explicit solution of the correlator equations (^), since then 
the correlator equations represent nonlinear coupled difference equations rather than a system of simple algebraic 
equations as in (H). The solution of the difference equations are uniquely determined by the boundary conditions 
induced by the two walls. Thus the pij will depend on the whole set of occupation probabilities {pi, . . . ,pm} in 
contrast to (||), where Pi^i+2m — fiiPi, ■ ■ ■ iPi+2m)- The best we can do for I > (2to + 1) is to reduce the problem to 
an equation for a single-site function hs = I — a^^^ /b^J^\ which is the probability that the site s is not occupied by a 
rod centre under the condition that the sites s — (/ — l),...,s — 1 are also not occupied by rod centres. Using /i^, one 
obtains the following set of equations equivalent to 



h, = 



l^r=s-{l-l) Pr + l^r=s~{l-l) l^k=r+2mPr,k 
l^r=s-{l-l) P»- + 2^r=s-(i-l) l^k=r+2mPr,k 



(i-/.,)e--^^nr 



, -1 



n 



i+l-l 
k=n 



{-^-h^)Y{ 



i+l-l 



^+Y:;^^+2^{i-K)e-<-Yv:i-'K^ 



(9) 
(10) 

(11) 



Solving ( |ll| ) for /i^, and inserting this solution into (10) gives the pi^j in terms of the ps. Inserting this solution then 
into (0) yields a closed system of equations for the ps ■ 

At this point we like to note that it is easy to generalise the method to lattice gases with internal degrees of 
freedom, where the occupation numbers can assume q different values x", a = 1, . . . , q (here and in the following the 
Greek letters have to be understood as superscripts). To be specific let us consider the (generalised) Potts model 



3|, where u"^ = u"''(5j,i+i. Following the procedure described above we find (with — '¥^[x1 



Q/3 



Ps~l,s — 



M f q q 



E 

13=1 

q 



Xlp"fi,,log^^ + log(l- 



Ps-l 



Z^a=l Ps-1 

Ps-1 



E' f a Q/3 \, /Ps ~ 2^f3=lPs-l,s\ 

".=1 p=l ^ 

q "1 

-, , aP \ , /-, Ps^ l^a,l3=lPs-l.s\ I 
l-Ps-l-Ps+ E Ps-l.s) l0g(l ) ' 



(12) 



and for the correlator equations. 



Q/3 

Ps-l.s 



1 -Ps 



Z^7,5=l Ps-1 



(13) 



The structure equations read 



EUp"-iJ{ps-E' 



S=l Ps,s+1 



12j,S=l Ps^l,s 



I-Ps-P. 



s+1 



■y.6-- 



--\Ps,s+\ j 



(14) 



independent of the position s these equations reduce to the ones previously derived by Percus |2 

for the pair-potential in eq. (H) and the associated eqs. 



For u^'' = u"^ 

Finally, we discuss the continuum limit of eq. (|^) for the pair-potential in eq. and the associated eqs. (|9|-|11D 
The continuous system of size L is composed of hard rods of length u that interact via a potential v{x, y) {y > x + a) 
exhibiting a finite range k < 2a. With p{x) being the density of rod centres and p{x, y) the two-point density, we 
obtain [Ea 




L pL — (T pK 

dyp{y)u{y)+ dy j dxv{y,y + x) p{y,y + x) 



dy\ 



JQ Ja 

P{y) - / dxp{x,y) 



log 



P{y) - Sy-l dxpjx, y) 



1 " /y-K dxp{x) + Jll_^ dx J^^^ dzp{x, z) 



-1 



+ dy dxp{x,y) log( 
Jo Jy-K 



The set of eqs. (p[]ll]) becomes 

h{y) 

p{x,y) 



p{x,y) 



Piy) - Iy~i^ dzp{z,y) 



-1 



(15) 



exp[-u{y)] = 



1 - ly-. p{z) + dx /J^^ dz p{x, z) 

p(x)h{y) exp[-'0(a;, y)\ exp (/J'^'' dz h{z) 
1 + J!^^^ dz h{z) exp[— {)(a;, z)] exp ^ /J^'" du h(u)j 

%)exp(/;+^dz/i(z)) 



1 + ly+a dz h{z) exp[-i)(a;, z)] exp(/J du h{u)) 



(16) 



(17) 



(18) 



Here h[y)dy is the probability that the interval (y, y + dy) is occupied by a rod centre under the condition that there 
is no rod centre in the interval (y — k, y). An explicit solution for h{y) and hence for p{x, y) can readily be obtained 
for the sticky-core fluid [p6[, and inserting this solution into ( p^ ) yields the functional derived in ||^. 

The functional in eq. (|5|) contains both the density and the density correlations, and to express the density 
correlations in terms of the density p{x) one has to solve a nonlinear integral equation for h{x). This form of the 
functional reflects the ad- hoc type of ansatz used in the weighted density functional formalism ||27| , which turned out 
to be one of the most successful methods in the theory of inhomogeneous fluids. The smoothing of the density on a 
local scale by a "weighting function" is in particular important to describe more subtle properties, as, for example, 
the wetting behaviour in the presence of confining walls . 

It would be desirable to establish a systematic procedure for deriving density functionals in higher dimensions 
d > 1. The method presented in this work can straightforwardly be generalised to higher dimensions as will be shown 
in forthcoming work. To keep the problem tractable, some approximations are needed in d > 1, but the resulting 
functionals will by construction satisfy the requirement of becoming exact upon dimensional reduction to d = 1 . In 
comparison to the methods applied in previous work |^] , the procedure based on the Markov chain has the advantage 
that the functionals are obtained directly and no functional integration of the structure equations is needed. 
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